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(N ! Abstract 



We study Chern-Simons number diffusion in a SU(2)-Higgs model with CP-odd dimension- 
eight operators. We find that the thermal average of the magnitude of the velocity of the 
Chern-Simons number depends on the direction of the velocity. This implies that the dis- 
tribution function of the Chern-Simons number will develop an asymmetry. It is argued 
that this asymmetry manifests itself through a linear growth of the expectation value of 
the third power of the Chern-Simons number. This linear behavior of the third power of a 
coordinate of a periodic direction is verified by a numerical solution of a one-dimensional 
Langevin equation. Further, we make some general remarks on thermal averages and on 
the possibility of the generation of the baryon asymmetry in a non-equilibrium situation 
due to asymmetric diffusion of the Chern-Simons number. 
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1 Introduction 
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An important cosmological observation is the matter-antimatter asymmetry in the present 
universe. This asymmetry may be quantified by the baryon-to-photon ratio, whose obser- 
vational value is Q 

— = (1.55 - 4.45) x 1(T 10 , (1) 

n 7 

with rig the baryon number density and n 7 the photon density. An interesting question is 
how this asymmetry was generated in the early universe. In his 1967 paper Sakharov p| 
notes that baryon number can only be generated when 

1) baryon number is not conserved, 

2) the transformations C and CP are not symmetries, 

3) there is a departure from equilibrium. 

In 1976 't Hooft || showed that in the elcctroweak theory the first requirement is 
satisfied. The non-trivial vacuum structure of the SU(2)-gauge theory in combination with 
the anomaly equation implies that a change in Chern-Simons number, NqSi is accompanied 
by a change in baryon number, B, 

B(t)-B(0) = 3[N CS (t) - JVos(0)] 

r.2 ft 



3^ 
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., , dt' / d 3 xF« v F^ a , (2) 
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with g the SU(2)-gauge-coupling, F£ v the field strength, and F^ va its dual. At zero tem- 
perature changes in the Chern-Simons number involve instanton processes. The rate of 
these processes is very small ||. However, at high temperature the rate is much larger. As 
has long been recognized this opens up the possibility that the baryon asymmetry was gen- 
erated at temperatures of about 100 GeV when the (elementary) particles that made 
up the plasma may be described by the electroweak theory. 

However, within the standard model the amount of CP-violation is probably too small to 
account for the observed asymmetry ([l]) || . One way to deal with this is to include effective 
non-renormalizable operators that parameterize the CP-violation of a more fundamental 
theory. The lowest-dimensional operator in a SU(2)-Higgs theory is $ (t>F^ v F^ va , with <j> 
the Higgs field. It may be included in the effective action as 



o[6] = 3 ff 2(5 cp f ,p X (tf6F a F^ a 
32tt 2 M 2 J 9 9 ^ ' 



(3) 



with M a mass and S^ P a coefficient that ideally may be derived from the fundamental 
theory. In the current investigation dimension-eight operators will play a more important 
role. They may be included in the action as 

S [8] = / d 4 x [6h P (D P <t>)HD P <t>) ~ S 2 CP F b pa F^ b + $ P (^W] F^F""*, (4) 

with Dp the covariant derivative. 

In this paper we study the combination of the first two of Sakharov's requirements, 
although we were unable to refrain from making some remarks on the inclusion of the 
third one in our study. That is, we study the motion of the Chern-Simons number in the 
SU(2)-Higgs action extended with (||) and @). It is well known that in a pure SU(2)-Higgs 
theory (without extra CP-odd operators and without fermions) the long-time behavior of 
the Chern-Simons number may be viewed as a diffusion process in one dimension, that can 
be specified by the expectation values 

(N cs (t) - N cs (0)) = 0, (5) 
([N CS (t) - iVcs(O)] 2 ) = 2VT sph t, (6) 

with V the volume and r sp h the sphaleron rate. The sphaleron rate plays an important 
role in scenarios for electroweak baryogenesis, since it determines the rate of the baryon- 
number violating processes. Much effort has gone into the determination of the sphaleron 
rate, both in the symmetric phase as well as in the broken phase |p| , [[lOf , |]ll|] . 

Here, only the broken-phase sphaleron-rate will be needed 

r sph - ktV^, (7) 

with the sphaleron energy -E sp h, and k a dimensionless coefficient that depends on the 
temperature, the Higgs mass and expectation value, and the gauge coupling. The physical 
picture that underlies (0) is that of classical transitions from one (classical) vacuum to the 
next over a potential barrier of height E sp \^. On the basis of this picture, the gauge-Higgs 
dynamics will be treated classically in this paper. 

The idea that the inclusion of CP-violating operators may have an interesting effect 
on the diffusion of the Chern-Simons number stems from the fact that the Chern-Simons 
number itself is CP-odd. Therefore the motion over the barrier towards positive Chern- 
Simons numbers may differ from the motion in the direction of negative Chern-Simons 
numbers. It is clear then, that the distribution function of the Chern-Simons number does 
not need to develop symmetrically. Or put differently: there is no symmetry argument 
why expectation values of odd powers of Ncs(t) — ^cs(O) should vanish. In fact, we will 
show that an asymmetry in the distribution function will develop indeed (as may have 
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been inferred from the title). That is, an initially symmetric distribution function will 
develop an asymmetry. In the presence of the CP-odd operators in (Q) an equation for the 
expectation value of the third power of Acs (t) — Acs (0) needs to be added to equations (^|) 
and @ to characterize the diffusion. We will argue that this expectation value increases 
linearly in time in the broken phase: 

({N cs (t) - N cs (0)f) ~t, (8) 

A more detailed conjecture is given in equation (fts|). 

We end the introduction with a short outline of the rest of the paper. In the next 
section the dynamics of the gauge-Higgs system will be projected on a one-dimensional path 
through the configuration space. This will simplify the analysis of the rest of the paper. In 
section || we will calculate the thermal average of the velocity when the system moves in the 
direction of positive or negative Chern-Simons numbers. In section ^ a conjecture for the 
time dependence of ([JVcs(i) — ^Vcs(O)] 3 ) is presented, where it is argued to grow linearly 
in time. This linear growth is also present in a simple random walk model as is shown 
in appendix Section ^| contains some general remarks on statistical averages. In the 
following section, we derive a simple stochastic equation that is useful to verify some of the 
analytic results. Also some numerical solutions are presented. In section ^ we discuss how 
asymmetric diffusion of the Chern-Simons number may yield a non-zero baryon number in 
a out-of-equilibrium situation. Also some comments are made on the constraints on such 
a model. We end with a summary in section pi 



2 Projection of the dynamics of the gauge and Higgs 
fields onto one special direction in configuration space 

To simplify the later study of the effect of the CP-violating operators in (|^) and (||) on the 
diffusion of the Chern-Simons number, the gauge-Higgs-field dynamics is projected on a 
single path through the configuration space. In order to include sphaleron transitions, we 
use a path that runs from a classical vacuum to the sphaleron configuration and further 
to the next vacuum. Manton [ fl2"| defined such a path in his proof that a non-contractible 
loop exists in configuration space of the SU(2)-Higgs model. This is the path that will 
be used in the following. The path of Manton is not the minimal-energy path, which 
was constructed in |Tq] . But the precise path will not be important for the following 
rather general arguments and we expect that the results will suffice as order of magnitude 
estimates. 

We parameterize the path of Manton by the time-dependent coordinate 9. The fields 
are given in the radial gauge (diA 10 - — 0) by 

gauge A^a a = ^f(p)[d^U(e)]U-\e), (9) 



Higgs 



4> = ±V2v lh(p)U(Q) fj)-i[l- h(p)} cos 6 ( J H , (10) 



with v the Higgs expectation value and p = gvr, where r = \x\. The O-dependent SU(2)- 
matrix is given by 

The functions / and h satisfy the boundary conditions 

/ ->■ h^O r^O, 

f->l r^oo. (12) 
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The form of the fields (|^), ( |l0| ) is a non-static generalization of the fields considered in 
(T^j . When O is chosen to be time- independent the fields given by ([)]) , (|l0[ ) can be identified 
(up to a global U(l) gauge transformation) with the fields (3.17) in ]12| (the coordinate 
corresponds to fi in Q). 

For the functions / and h the Ansatz b of Klinkhamer and Manton |Q will be used 



m 

Hp) 



A(A+4) 

1- ' exp[i(A-p)] 



A+4 

oB + 1 
B(aB+2) 



p<A 
p> A 



P 



1 



B 



1 

a-B+2 p 



exp 



\{B-P)] 



P<B 
P>B 



(13) 
(14) 



with a = (X/2g 2 )^ , where A is the Higgs self-coupling constant. The parameters A, B were 
determined by minimizing the energy for the static field configuration at = it/2 in fl4|| , 
in order to obtain an estimate for the sphaleron energy. In this way the parameters depend 
only on a. We take a = 1, for which the parameters have the values A = 1.15 and B = 1.25 
& 

Now that the dynamics has been restricted to the path described by (g) and (|l0|) we 
may rewrite the SU(2)-Higgs action, S, in terms of the coordinate 6 



S[Q, Q] = ^j- I ,11 



(ax + a,2 sin 2 0) 



e 2 



M 2 



(03 sin 2 6 + a4 sin 4 9) 



The CP-odd action (^) in terms of reads 

47TU 2 



@M = ^r J dt [Mcp + fe^cp + (Mcp + &4<5cp) sin 2 6] 6 3 sin 2 9, 



(15) 



(16) 



where total time-derivatives have been neglected. 

The action (||), that includes the dimension-six operator tfttpFF, gives only a 
total time-derivative, just as the operator 4>) 2 F^F* 11 " 1 in (^). Therefore the inclusion 
of will not affect the motion along the path parameterized above. When the full 
field-dynamics is considered, the dimension-six operator may have an effect on the motion 
along the Chern-Simons number direction in a non-equilibrium situation, see e.g. fl5|| . 
In equilibrium however, the system oscillates around one vacuum or, during a transition, 
around the minimal-energy path that connects two classical vacua (the path that we have 
approximated by (JsJ) and (|lC|)). Since, the oscillations average out in equilibrium, we do 
not expect the inclusion of^ 6 ' to affect the diffusion of the Chern-Simons number. 

The coefficients a±, a,2, 03, 04, &i, 62, b 3 64 are given by the integrals 



"2 



«3 



0,4 



dpp 



2.03, 



dp 



dp 



dp 



8f(l-f) 2 + -(l-h 2 )(l-f)p 2 

^d p ff+ l - P 2 {d p hf + h 2 {l-ff 
-2fh(l - /)(1 - h) +f(l-h) 2 ] =1.41, 
8f 2 (l-f) 2 /p 2 + 2fh(l-f)(l-h) 

-f 2 {l-h) 2 + ^p 2 (l-h 2 ) 2 =0.70. 



= 2.53, 



(17) 
(18) 

(19) 

(20) 
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h = 2 / ^(^-/) 2 (<9 P /)/(l-/) = 0.09, (21) 



b 2 = 9 / dp{d p fYf{l- f) = 0.10, (22) 



6,3 = | / '//'CI ■/riUW'l/d - /id -./^ =n.l<). (_>:-!) 



f°° 1 

6 4 = 72 / dp-(d p /)/ 3 (l- /) 3 = 0.16. (24) 
Jo P 

Let us shortly discuss the actions ( |l5| ) and ([l6]). In the action (|l5|) we recognize the 
potential barrier fll^] , 

V(e) = -y^ [03 sin 2 + a 4 sin 4 9] , (25) 
between different vacua. The sphaleron energy reads 

1 4:7TV Attv 

E sph = V(Q = -tt) = a 3 + a 4 = 2.11 . (26) 

2 5 5 

The sphaleron energy determines the exponential suppression of the sphaleron rate (Q). 
Another quantity that enters the sphaleron rate is (the real part of) the frequency of the 
negative mode, w_, at the sphaleron configuration ||. From the action ( |l5|) it may be 
determined as 



w_ = 



Oi + a 2 



« 0.61( 5 «) 2 . (27) 



In agreement with the order of magnitude estimate uj- ~ in Q . 

The action (|l6|) is odd in 0. This can be understood by realizing that a CP-transformation 
in terms of O is 

e^-e, (28) 



since the action (Jl|) is CP-odd, so should (|16[). 

Now that the SU(2)-Higgs action including the effective CP-odd operators has been im- 
mensely simplified, namely to the one-dimensional action ( |T5| ) + (|l6|), the question: "does 
the Chern-Simons number diffuse asymmetrically under the influence of the dimension- 
eight operators in the action (0)?" may be cast in the form: "does the factor 3 in ( [l6|) 
lead to an asymmetric diffusion of 0?". What is meant by "diffusion of 0" may require 
some explanation. From equations (P) , (|i"o|) , and (pr|), one sees that = and = 2-7T 
correspond to the same (static) gauge and Higgs field configuration. Even = and 
= 7r correspond to the same physical state, since they are related by a simple gauge 
transformation. Hence, one may view configuration space as a circle with circumference tt. 

This is the same for the full gauge-Higgs fields. When a transition over the sphaleron 
barrier is considered, the system starts and ends in the same state (since the different vacua 
are related by a large gauge transformation). Nevertheless the winding number is changed: 
ANqs = 1. The (change in the) winding number is the physical quantity of interest, since 
it is related to the (change in the) baryon number through (||). The relation between a 
change in the Chern-Simons number and is (see Appendix A) 

A^cs (t) ~ A cs (0) = - / dt'Q sin 2 0. (29) 

7T 







Instead of considering the system on a circle and keeping track of the winding number, 
it is convenient to "unwind" the system and consider the system on an infinite line £ 
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Figure 1: The distribution function of O in a periodic potential. On the left it is viewed on the 
circle 6 [0, 7r[, then the distribution function is static. On the right it is unwound to the infinite 
line G] — co, oo[ accounting thus for the winding number naturally. On the infinite line the system 
is dynamical and the distribution function diffuses. 



] — oo,oo[ and follow the dynamics of 9 (see figure 0). Thus, we consider the diffusion of 
6] — oo, oo [, where we keep in mind that actually we consider the winding number on a 
circle. 

To discuss the diffusion of in the next section it is convenient to derive the Hamilto- 



nian. Firstly, we write the Langrangian corresponding to (|15|) and (16) as 



L = z 1 (e)e 2 + z 2 (e)e 3 -^(6), (30) 

with 

47T 

zi(0) = -«-(oi + a 2 sin 2 0), (31) 

z 2 {@) = 1?r [6 1 ^p+McP + fe^p+Mcp)sin 2 e]sin 2 e. (32) 

Clearly the kinetic energy in ([50|) is unbounded. And, for instance, a (CP-even) 4 is 
required to put a lower bound on the kinetic energy. Such a term should come from 
a higher (non-renormalizable) operator, so it is expected to be suppressed compared to 
the other CP-even 2 -term. Therefore the details of the 4 -term are unimportant for the 
following calculations. In fact strictly to first order in 8q P and (which are assumed 
to be small), one can work with the unbounded Lagrangian (|3C|) (and the unbounded 
Hamiltonian (^J)), although it may be kept in mind that a 4 -term (or higher even powers 
of 0) are required for a bounded energy. We assume that these higher-order CP-even terms 
will make the Lagrangian convex, so that we can go to a Hamiltonian description. The 
conjugate momentum then reads 

p = 2^(0)0 + 3z 2 (0)0 2 , (33) 

and the Hamiltonian 

This is the Hamiltonian that will be used for the calculation of statistical averages in the 
next section. 

The point of this section was to show that the CP-odd operators in (|4|) projected 
onto the one-dimensional path in phase space gives a p 3 -term in the Hamiltonian and 
a 3 -term in the Lagrangian, and to give reasonable estimates for the coefficients in the 
projected action and Hamiltonian. The main subject of this paper, the effect of the CP-odd 
terms on the dynamics, will be discussed in the next and following sections. 



3 Velocity expectation values 

We study the effect of the action (jij) on the diffusion of the Chern-Simons number. For 
this the Hamiltonian ( |3~4| ) will be used. We will work to first order in the coefficients Sq P 
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and (5q P (and hence, to first order in z 2 ). So the question is what does the diffusion of © 
look like to first order in <5q P and <5q P . In particular, will the distribution function of 
develop an asymmetry? 

The first thing one may notice is that (the magnitude of) the velocity differs for positive 
or negative 0. More precisely, at a given momentum (p) or at a given energy (H) the 
magnitude of the velocity depends on its direction. This implies that the motion over the 
barrier towards negative Chcrn-Simons number differs from the motion towards positive 
Chern-Simons number. Hence, an asymmetry is expected to develop. 

It is useful to define a quantity that "measures" the asymmetry. A simple and useful 
quantity is the mean velocity, at a certain position, when moves either to the right or to 
the left. We define 

v^e) = (\Q\H(e)6(e-e))/(H(Q)5(e-e)), (35) 
v l (e) = (\@\n(-&)s(@-e))/(n(-e)s(e-e)), (36) 

where the brackets denote a thermal average, and Tt the Heaviside function. When the 
difference 

Av(0) = v T (0) -V(0) (37) 

is non-zero, an asymmetry will develop. Consider, for instance, the evolution of the distri- 
bution function of 0. When the particles^] that move to the left do it faster (slower) than 
the particles that move to the right, the tail of the distribution function to the left extends 
more (less) than to the right. Hence, the distribution function develops an asymmetry. In 
particular, we expect expectation values of odd powers of to become negative (positive) 
when the distribution function is initially symmetric and thermal. 

In the calculation of the velocities ( |35| ) and (^6j), firstly one may remark that their 
numerators are equal. This follows from the fact that the flux through a point 

(05(0 - 0)) = (|e|W(0)*(e - 0)) - (\@\H(-Q)5(Q - 0)) (38) 

vanishes. This may be seen from 

(05(0-0)) = (d p HS(@ - 0)) 



poo f K /2 

= Z- 1 / dp d<dd p H5{e - 0)e- /3ff 

Joo J-n/2 



I-tt/2 
poo /' 7r /2 

= -Z~ X T \ dp dQ6{e - Q)d p e-P H = 0. (39) 

Joo J-n/2 

Here we used that exp— /3H vanishes at p = ±oo, for which we had to keep in mind that 
a p 4 -term in the Hamiltonian is required to make it bounded and convex. (Alternatively, 
strictly working to first order in <5q P and S^p, the momentum integrations in ( |39"| ) also 
vanish.) The normalization factor is defined as usual, see (|5l]). It will drop out of the 
calculation of the velocities and v+. The fact that the flux vanishes, implies that the 
expectation value of remains constant: 

(0(t) - 0(0)) - 0, (40) 

as will be discussed more fully in the section S. Hence, asymmetric diffusion will affect 
only expectation values of higher odd powers of 0. 



It is convenient to think of the distribution function as a normalized sum of a lot of single particle positions. 
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The result for the numerator on the r.h.s. of ( |35[ ) (equal to the numerator on the r.h.s. 
of ph) is 



/oo fx f2 

dp / d<dd p HH( P )5(@ - 9)e- /3ff 
-oo J-rr/2 



1^-/3^(0) 



Z-'Te 



(41) 



The denominator of (35) is given by 

(w(e)*(e-6)) = z^-w / d pe ^^rw) p + ^w p 

,2a(e) _ 



7r2i(e)T + T 



to first order in z%. The denominator of fl3q ) yields 

(W(-0)<5(0 - 0)) = Z- 1 \J irzi(0)T - T 
Hence, up to first order in 5q P and Sq P we have 



z 2 m 



,1(1) 



(0) 



T 



7^(0) 



and 



(+: 



2T 



T 



7r[*i(0)] : 



r^(6), 



Tr^(e)]^ 



(42) 



(43) 



(44) 



(45) 



The relative velocity difference, 



A^ rel (9) 



Av(0) 



v1(S) 



-v1(Q) 



[zi(B)} 3/2 



(46) 



gives a (dimensionless) measure for the effect of CP-violation on diffusion. 

In conclusion, we have established that, when 0^0, the magnitude of the thermal 
expectation value of the velocity depends on the direction in which moves. 

Let us consider what the difference in thermal expectations values of the velocity implies 
for the evolution of the distribution function of and p, P(p, 0, t). When the momentum, 
p, is thermally distributed independent of the position 0, the velocity difference is present 
in the entire (0-)space. Now it is not hard to imagine that, when the particles move faster 
to the left than to the right, the tail of an initially symmetric thermal distribution function 
will (start to) extend more to the left than to the right. Therefore, a symmetric thermal 
distribution function will not remain symmetric. The diffusion develops asymmetrically 
when CP-odd terms like 3 are included in the lagrangianQ. When we translate this result 
back to the gauge-Higgs system it implies that the inclusion of the CP-odd operators in 
(JiJ) leads to an asymmetric diffusion of the Chern-Simons number. 



2 Strictly, the argument given only implies that P(Q, p, t) — P(— O, p, t) cannot remain to hold as time evolves. 
In principle, it might still be possible that f dpP(<d,p,t) — f dpP(—<d,p,t) remains true. Since, in the broken 
phase at least, we expect the momenta to be thermally distributed independent of the position, O, we do not 
consider this possibility further. 
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4 Asymmetric diffusion: a conjecture for large times 



In the previous section, it was shown that there is an asymmetry in the average velocity ( |46| ) , 
which implies that the distribution function will develop an asymmetry. Hence, for instance, 
the expectation value of the third power of the Chern-Simons number, ([Ncs(t) — Acs(0)] 3 ), 
will be nonzero. In this section, we will argue that this expectation value will grow linearly 
in time in the broken phase of the SU(2)-Higgs model. 

We assume that after a transition over the sphaleron barrier the system thermalizes 
before a following transition. This implies that the transitions are uncorrelated. This 
assumption is reasonable when the height of the barrier, i? sp h, is large compared to the 
temperature, T, and the motion over the barrier is damped sufficiently. (In our description 
in terms of 0, this would require O to be not-too- weakly coupled to the other modes in 
the plasma.) 

On the basis of this assumption we can establish the time-dependence of the third 
power of the Chern-Simons number. We have established in the previous section that an 
asymmetry will develop in the distribution function starting from an symmetric thermal 
initial distribution function. Hence, after a time At the expectation value ([iVcs(At) — 
iVcg(O)] 3 ) has a nonzero value. During the time At the distribution function has spread over 
different vacua. Since we assume that the distribution is thermal in each vacuum, during 
the time from At to 2At, from each vacuum the distribution will diffuse asymmetrically in 
the same way as it did in the time interval from time to time At. The thermalization in 
each vacuum implies that in each time interval At, from each vacuum the same diffusion 
process takes place (relative to that vacuum). From this argument follows the result that 
the expectation value of the third power of the Chern-Simons number in the presence of 
CP-odd dimension-eight operators in the broken phase of the SU(2)-Higgs model grows 
linearly in time. That the repetition of the same asymmetry in the diffusion process in 
each time interval leads to a linear growth, may be seen, for instance, from a random walk 
model, as is shown in Appendix B. 



In section 5T, we show that in general (without the above asumption) the expectation 
value of a third power of a coordinate x, (x 3 ), either stays constant or grows linearly in 
time in thermal equilibrium. A stochastic model that will be introduced later, shows that, 
with potential barriers, (a: 3 ) grows linearly in time, and that without barriers (x 3 ) goes to a 
constant. This supports the view that barriers are required to ensure that the asymmetry 
in the diffusion process is independent of the position (which vacuum the system is in) 
and time, from which the linear behavior follows. In the symmetric phase, different vacua 
are not well separated by a barrier, therefore it is for us not possible to determine the 
time-dependence of asymmetric Chern-Simons number diffusion. 

To obtain a more detailed conjecture for diffusion in the broken phase, we extend our 
one-dimensional model by coupling to the other modes (that form the heat bath at 
temperature T) as in the full SU(2)-Higgs model. In this way the rate for transitions over 
the barrier is given by (<7t>) -3 r sp h, with (gv)~ 3 the volume of the sphaleron. Then, we 
expect 

(0 3 (t)) ~ (^)- 3 r sph tA WreliSph , (47) 

with the relative velocity difference ( pi] ) at the sphaleron configuration Aii rc i,sph = Aw ro i(7r/2) 

By viewing the entire space as made up of blocks of volume (gv)~ 3 in each of which 
a 0-coordinate is diffusing, one can translate ( f47| ) into an expectation value for the third 
power of the Chern-Simons number 

([A cs (t) - A^ C s(0)] 3 ) = C yr sph tAtw, sph , (48) 

with V the volume and c a constant of order one. In going from (]47|) to ( ftsj ) we have in 
a very non-sophisticated manner included the zero mode for translation invariance. (It is 
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conventional to include rotational zero modes in T sp h and exclude the translational zero 
modes.) A better way would have been to include the zero modes in the parameterization 
@ and (0). However our aim is to investigate the occurrence of asymmetries for which 
the zero modes play no essential role. 



5 Remarks on statistical averages 

The Chern-Simons number diffusion is a non-equilibrium process. Still we, as many others, 
see e.g. 0], use a thermal (equilibrium) average to calculate the mean or average velocity 
of the Chern-Simons number (of actually). In this subsection, we will briefly consider 
the use of statistical averages. In the next subsection, we return to dynamical issues and 
consider the possible time-dependence that an expectation value can have. 

Consider a particle in one dimension coupled to a heat bath. In general one expects 
that a probability distribution function, P(p, x, t), of the position, x, and momentum, p, 
in the long time limit, goes to a thermal distribution function: 

P(p,x,t) ^-i e -/3H(p.*) j (49) 

with Z~ x the normalization factor. Then, the long-time limit of expectation values can be 
calculated using the thermal distribution function. However, in the case that the Hamil- 
tonian is periodic, the thermal distribution function is not normalizable. A well-known 
consequence is that (x 2 (t)} in a flat or periodic potential does not have a thermal limit- 
value at large times (since the thermal average is not defined) . Instead it grows linearly in 
time. Also thermal averages of other positive powers of x cannot be calculated. In partic- 
ular, since (x n ) is not well-defined in equilibrium, one cannot conclude from a symmetric 
potential, V(x) = V(—x), that for odd n it should vanish. 

Nevertheless there is still a lot that may be calculated using a thermal distribution. This 
is based on the notion that the long-time limit of the distribution function in a periodic 
potential with period n satisfies 

*£P(p,x + mr,t) ^ Z^e-f B <*"\ (50) 

71 

where Z~ x is the normalization factor, with 

/7T /2 P OO 

dx / dpe~ (iH . (51) 
-7T/2 J-oo 

Equation ( |50| ) expresses the fact that when the points x + mr are identified, that is when we 
view the dynamics of the system on a circle with circumference it, the system thermalizes 
as usual. There are two different ways to look at a periodic system. When the system 
is considered in terms of the coordinate x G] — oo, oo[, the system is dynamical. There is 
diffusion, transitions to other classical vacua, etc. Whereas, when we view the system as 
living on a circle, with coordinate <p = x mod7r, the system is static after thermalization 
(of course, the dynamics returns when the winding number is considered). 

Thermal averages can be calculated of functions, f(p,x), that are either independent 
of x or 7r-periodic in x. In the long-time limit these functions have a thermal limit value: 

/■OO p^TT />OG 

dx / dp f(p, x)P{p, x, t) 1 Z- 1 \ dx dp f(p, x)e-P H <*>*) . (52) 



In particular the average velocity can be calculated with a thermal distribution (as we did 
in section ||). Since the average velocity vanishes 

(x) = 0, (53) 
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one expects that (x(t)} is constant (in the long-time limit). There exists no such argument 
for expectation values of other odd powers of x. For instance, the expectation value of the 
time-derivative of a; 3 , (3ia; 2 ), is itself not well defined. Hence, (x 3 ) does not need to vanish 
and has no thermal limit-value. Indeed, as we have argued earlier, it may grow in time. 

5.1 Possible power laws 

Next, we will determine on general grounds the possible time-dependence of (x 3 (t)) in a 
periodic potential in the long-time limit. We assume that it behaves as a power law 

(x 3 (t)) =C[I\ a ]1?, (54) 

with C and p constants. The brackets denote a classical average, that is, an average over 
initial conditions weighted by an initial distribution function P\ n . The constant C may 
depend on the initial distribution. 
We will use that 

(x(t) - x(0)) = 0, (55) 

see (^). It is convenient to choose the origin such that (x) vanishes at the initial time. 
Since then (x(t)) = for all times, and the disconnected contributions to (x 3 (t)) vanish. 
To determine the possible values for p it is useful to consider 

([x(t 2 ) - x(h)] 3 ) = (x 3 (t 2 )) - 3{x 2 (t 2 )x(h)) + 3{x(t 2 )x 2 (h)) - (x 3 (h)). (56) 

When t 2 — t\ — > oo, we can apply ([si]) to the left hand side 

([x(t 2 ) - xih)} 3 ) = C[P(h)}(t 2 - hf. (57) 

The constant C[P(t±)] depends on the distribution function at time t\ (instead of the initial 
distribution function). A small subtlety is that P{t\) should be considered as a distribution 
function of the shifted coordinate y(t) = x(t) — x[ti) instead of x{t). 

Consider the correlation functions (x(t 2 )x 2 (ti)) and (a; 2 (t2)^(^i)) on the right hand 
side of (|5^) in the same limit t 2 — t\ ^> oo. When the system thermalizes it follows from 
the fact that the expectation of x goes to a constant (|5j]) that the correlation functions 
([x(t 2 ) — x(ti)]x 2 (ti)) and ([x(t 2 ) — x(ti)] 2 x(ti)) go to a (small) constant in the limit 
t 2 — t\ — > oo. Therefore, when p > 0, we have in this limit 

(x(t 2 )x 2 {h)) -» (x 3 ^)), (58) 
(x 2 (t 2 )x(h)) -> (x 3 ^)). (59) 

When both t 2 and t\ are send to infinity also, ([561 ) gives 

C[P(h -» oo)](t 2 - hf = C[P m ]t p 2 - C[P m ]t\. (60) 

This equation has two possible solutions 

p = 1, with C[P{h oo)] = C[P m ] = C, (61) 
p = 0. (62) 

Hence the expectation value (x 3 (t)) either grows linearly in time or stays constant. When 
the expectation value (x 3 (t)} goes to a constant in the long-time limit, it is not so that 
this constant should be zero. The constant can be non-zero when it depends on the initial 
distribution. (In appendix C we consider a stochastic equation where (x 3 ) indeed goes to 
a non-zero constant in the long-time limit.) 
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It is reasonable to expect that a similar reasoning can be applied to connected correlation 
function of higher powers of x, with the result that these too are either constant or linearly 
growing in the long-time limit. Then, in the case that (x 3 (t)} grows linearly in time, the 
expectation values of higher powers of x are dominated by their disconnected parts: 

<* 2 " +1 (*)> « 2 n%l^ 1 y ^(t))[(AtW- 1 - (63) 

For a linearly growing expectation value of a; 3 , we get (x 2n+1 ) ~ t n . 

To summarize, on general grounds it has been shown that the expectation value of 
(x 3 (t)) goes to a constant or grows linearly in the long-time limit. As we have argued in 
section ^ we expect that, in the broken phase, the expectation value of the third power of 
the Chern-Simons number grows linearly in time. It is encouraging that the result of the 
argumentation in section [| is consistent with the more general analysis presented in this 
section. Nevertheless, it is desirable to have a more realistic model that would enable one 
to test the conjecture that (x 3 ) in a periodic potential with high barriers grows linearly in 
time. In the next section we present such a model, namely a simple stochastic equation. 



6 A stochastic equation 

In section ^| a one-dimensional Lagrangian was derived for the motion of the Chern-Simons 
number. In this section we want to (re)introduce the effect of the (infinite number of) 
other modes [ [To) . They provide a heat bath at temperature T. The simplest way to mimic 
their effect is by the introduction of a damping term and a stochastic noise in the equations 
of motion. In fact, the motion of the Chern-Simons number in the symmetric phase is to 
leading order indeed determined by a stochastic (field) equation H. In the broken phase, 
the case that is of interest to us, a local damping term and stochastic force is probably only 
suitable for illustrative purposes, and not of direct quantitative interest for Chern-Simons 
number diffusion. Nevertheless the stochastic equation derived below will provide a simple 
and realistic model for asymmetric diffusion. 
Consider the Lagrangian 

L = \{x 2 + 8x 3 + dd 2 x 4 ) -V(x), (64) 

with S a small parameter and d a dimensionless coefficient. The x 3 is the time-reversal 
non- invariant term similar (but simpler) as the O 3 term in (^0|). The i 4 has been added 
so that the kinetic energy is bounded from below. We will demand that the Lagrangian is 
convex, which requires d > 9/24; we choose for the following d = 1. 
The equation of motion is 

x (l + 35x + 6S 2 x 2 ) = -d x V. (65) 
We now introduce damping and a stochastic force in the following way 

x(l + 36x + 66 2 x 2 ) =-d x V-jx + £, (66) 
with 7 the damping coefficient and £ a Gaussian white noise 

(£(*)) = 0, (67) 
WW)) = 2 7 TS(t-t'). (68) 

The introduction of the damping term (72;) looks standard. However, there is one sub- 
tlety. This may be made clear by going to the Hamiltonian formulation. The Hamiltonian 
corresponding to ( |64| ) (with d = 1) is 

H = i (p 2 - S P 3 + 8 2 p i ) + V(x) + 0(S 3 ). (69) 



12 



The Hamilton equations up to order S 2 including noise and damping in the same way as 
in (|66|) read 



x = v p =p-^Sp 2 + 26 2 p 3 , (70) 



3 

p = -d x V--fv p + £. (71) 

However the introduction the damping term as "fv p instead of jp is arbitrary at this point, 
and requires an explanation. 

One argument for the introduction of the damping term as "fv p goes as follows. In a 
microscopic derivation (for example using influence- functional techniques []l9|| ), one would 
find that integrating out the modes of the heatbath yields a memory kernel of the form 
J dt'H(t — t')x(t'). In a short-time expansion, this memory kernel will reduce to jx. This 
reduction is independent of the (possible complicated) structure of the kinetic term in 
the Lagrangian or Hamiltonian. Therefore we expect that in time-reversal non-invariant 
theories with complicated kinetic terms, a derivation from first principles would yield a 
damping term of the form 72; as introduced in ( |66| ) and (|70|) , (f7l|) . 

Perhaps a more convincing argument is given by considering the equilibrium distribution 
function. The equation for the time evolution of distribution functions, the Fokker-Planck 
equation, may be derived in a standard manner, as for instance in ]l6|] . It reads 

t P = d p frTdp + (d x V) + jv p ] P - v p d x P, (72) 

with P = P(p,x,t) the time-dependent distribution function. The static solution to this 
Fokker-Planck equation is found to be exp —f3H. This shows that damping is correctly 
introduced in ( p6| ) and (|7l|). 

The remarks of section |^ apply to the stochastic system ( |70| ) and (71). Namely, for 



periodic functions or functions that are only dependent on the momentum p, the thermal 
average determines the long-time limit (p2[)P|. For instance odd powers of the velocity have 
the following limit 

with n = 0, 1, 2, .... Similarly, it may be concluded that, for example, the long-time limit of 
odd powers of sin(2x) vanishes when the potential is symmetric and periodic with period 

7T. 

The stochastic equation ( |66[ ) introduced here will be used in appendix ^ for some sample 
calculations with potentials V = and V = \lJ 2 x 2 . 



6.1 Numerical simulations 

Perhaps the most relevant case for the stochastic model described above is when the poten- 
tial V is a periodic function similar to ( p5| ) describing a potential barrier. For simplicity we 
take V to be V(x) — V h sin 2 (a;) where V h is then the potential barrier height. The solution 
to this case is provided by numerical analysis of equations ( |70| ) and (|7l]) for a considerable 
amount of realizations. Average quantities are then obtained by averaging over all those 
realizations. The results from a numerical simulation are displayed in figure^ for the evo- 
lution of the system up to large times (t ~ 2500). One immediately sees from figure ^ that 
(x 3 (t)) grows linearly in time. In section |5| we obtained two possible long-time behaviours 
for (a; 3 ), namely (x 3 ) = constant or (x 3 ) ~ t. Therefore the numerical results support 
our arguments in favour of the latter for a periodic potential. The fact that it grows to a 

3 Also for periodic functions or functions that only depend on the momentum p the ergodic theorem implies 
that the statistical average equals the time average. Even when the time-reversal symmetry is broken jl7j|. 
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time (t) time (t) 

Figure 1: Evolution of {x 3 (t)) (left) and (x 2 (t)) (right) for S = 0.1, 7 = 0.3 and Vh = 4 (in units 
where T — 1) with 8 • 10 J realizations. Statistical errors of one standard deviation are also included, 
together with their time evolution (small pictures). These grow as t~ for (x n (t)). 



positive value is somewhat surprising since our equilibrium considerations showed that the 
velocity in the negative direction v^(Q) is larger than the velocity in the positive direction 



u T (9) for S C p > 0, so Av 



= IV — 7)1 



< 0. This behaviour for the asymmetry Av in the 



velocities is also found in our numerical simulations. It would seem to follow that (x 3 (t)} 
evolves in time to a negative value. This is also supported by the random walk model 
considerations of appendix |b].[|] The numerical results however indicate the opposite, which 
means that probably non-equilibrium effects determine the sign of growth of (x 3 (t)). It 
might also be that the sign depends on the details of the interactions, so that a different 
model could lead to a different sign. 

Expectation values of higher powers of x are dominated by their disconnected parts in 
the long time limit as argued in section [5]. Therefore, according to equation (|6^) one has 



/»(*) = 



„2n+l 



(*)> 



(2n+l)! 
T=l 



(74) 



We have explicitly checked this for (x 5 (t)) and (x 7 (t)} (figure (^) below). 





time (t) 

Figure 3: Evolution of f 5 (t) = 10( J^]^^ (left) and f 7 (t) 



time (t) 

i05<J(t)>*L^(t)>2 ( n 9ht). 



In the case of (x 3 (t)), its disconnected part (x 2 (t)} ■ (x(t)) is negligible. This is because 
(x(t)) is constant in time and thus equal to (x(0)) which was chosen to be (x(0)) = in 
our simulations. Actually in the numerical results |(x)| is of the order of 0.02 ± 0.1. The 



Nevertheless it is important to say that one cannot make a direct correspondence between the A in 
random walk model of appendix |b| and the 5 considered here in this stochastic model, not even for the sign 
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large-time limit of (x 3 (t)) in figure || corresponds therefore to its connected part. 

We have also checked in the simulations that velocity expectation values (v 2n+1 ) have 
the limit given by ( |73| ) . Thus we have verified that our stochastic model thermalizes in the 
way as expected in section ^J. 

7 On baryon-number generation from asymmetric Chern- 
Simons number diffusion 

A natural question to ask is whether the asymmetric diffusion of the Chern-Simons number 
may yield a non-zero baryon number in some out-of-equilibrium situation. This section 
contains some remarks on this question. In particular, we will present a simple way in 
which asymmetric diffusion of a particle, with coordinate x, will lead to a (temporary) 
non-zero expectation value of x itself. It is then argued that such a situation occurred 
during the electroweak phase-transition yielding a non-zero Chern-Simons number (and 
baryon number). In any electroweak baryogenesis scenario the two important questions 
are: Is the resulting baryon asymmetry sufficient to account for the observed asymmetry 
in the universe (|l|), and what are the requirements to prevent a wash out of the generated 
asymmetry? We will present some estimates to answer these questions. 

Let us reconsider the diffusion of a particle in one dimension with coordinate a:, that 
starts at x = 0. As time evolves the distribution function will spread. In the long-time limit 
the distribution will become Gaussian with on top a small asymmetry (assuming that the 
symmetry-breaking terms are small). We have argued that (x 3 (t)} grows linearly in time. 
To be definite, let us say it grows in the negative direction. Then the tail of the distribution 
in the negative- a; direction (-tail) will be larger than the tail of the distribution function 
in the positive- a; direction (+tail). This difference in the tails accounts for the negative 
values of (x 3 ) (and of expectation values of higher powers of x). That the expectation 
value of x remains zero is due to an asymmetry in the distribution function closer to x — 0. 
(A nice way to imagine this that the peak of the distribution function (the most probable 
position) moves in the positive- a; direction. The motion of the peak would give a growing 
(a;) except that the contribution of the tails is precisely opposite. For higher powers of x 
the tails dominate and a negative and growing value is the result.) With this picture of 
the evolution of the distribution function it is not hard to construct a system for which 
(a;) itself becomes (temporarily) non-zero. All that one has to do is to prevent the -tail 
to compensate for the positive contribution of the peak moving in the positive- a: direction 
(if the peak does not move in the positive a;-dircction there is still a positive contribution 
from the region around x — and the argument goes through unchanged). Consider the 
diffusion of this particle in a box with symmetrically placed walls at x = ±a. The -tail 
will hit the wall earlier than the +tail. Then the -tail can no longer compensate for the 
asymmetry of the distribution function close to x = (for instance the motion of the peak 
of the distribution function in the positive- a; direction). Hence, (a;) will grow and become 
non-zero. Eventually, when the system goes to equilibrium, (x) relaxes back to zero. 

Instead of putting the system in a box we could also have let the system evolve in a 
harmonic potential V — \ui 2 x 2 (to be superimposed on the periodic potential in which it 
diffuses). Basically by the same arguments as above it follows that the expectation value 
of x becomes temporarily non-zero. A slightly more general case that may be considered, 
is a system that is initially in thermal equilibrium in the presence of a harmonic potential 
Vi — \u\x 2 which changes at the initial time to V2 = \ijj\x 2 , with uj 2 > In the 
evolution towards the new equilibrium state, the expectation value (a;) will temporarily 
become non-zero. In the next section we will show in a short-time expansion that (a;) will 
become non-zero indeed, after such a change in a harmonic potential. 
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7.1 Short-time expansion 

We show in the following that an instantaneous change from a potential V\ = V pcl - + \io\x 1 
to V 2 = V-pc-r + \10\x 1 ', with V pcl a periodic potential, will lead to a nonzero expectation 
value (x). We will use the Fokker-Planck equation (f72|), with V 2 as the potential, for the 
time evolution of the distribution function P(p,x,t). The initial distribution function at 
time t = t- m is given by 

P ia (p, X ) = ^-l e -/3[i(f 3 -^ 3 +«V)+^W] ) (?5) 

and the expectation value (x) by 

(x(t)) = / dx / dpxP(p,x,t). (76) 



We calculate this expectation value in a short-time expansion 

^ = E -* x »r, (77) 



with coefficients 



Define the operator 



x n = jdx J dpx —P(p,x,t) 



(78) 



t=t h 



= d p frTd p + jv p + (d x V)]-v p d x , (79) 

such that the Fokker-Planck equation (|7|) may be written as dtP = OP, with potential 
V = V2. The coefficients x n may be calculated by 

x n = j ' dx J dpxO n P in (p,x). (80) 

The first three terms in the expansion ( f77| ) vanish 

x = {x) in =0, (81) 

xi = (« p > ln = 0, (82) 

x 2 = ((l-36 P )d x V 2 - 1 v p ~36 1 (T^p 2 )) in ^0, (83) 



where 



For X3 we find 



J dx J dp..P in (p,x). (84) 



a; 3 = 38(Ta£V 2 - (d x V 2 f)- m - (85) 

When V2 = V\ this gives X3 = 0. To first order in cj^ — w| the result for X3 is independent 
of the periodic potential, V^er, 

x 3 = 3S{luI - luI)T. (86) 

This shows that a change in the potential will lead to a (temporary) non-zero value for (x) , 
even if the initial and final potential is symmetric. 
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7.2 Numerical results 



To obtain the behaviour of the evolution of (x) for longer times we numerically solve the 
stochastic model discussed in section o, with the instantaneous potential change Vbcforc = 
Vbsin 2 (a;) + \to\x 2 — > t4ft e r = Vb sin (x) + \lo\x 2 with us\ and (J 2 chosen to be equal to 
0.1 and 0.01 respectively. The system is first let to be equilibrated in 14>efore and then the 
instantaneous potential quench is performed. The numerical results show that indeed a 
temporary non-zero value for (x) is obtained (see figure ^) . 




Figure 4: Evolution of {x(i)) (left) and (x 3 (t)) (right) after the instantaneous potential change 
Vbefore — ► V a ftcr (at t ~ 1000,) with statistical errors included. 



7.3 Baryon-number generation 

We now turn to the problem of baryon-number generation. The back reaction of the 
generated baryons on the Chern-Simons number may be described by the effective potential 
given by the free energy at a given baryon-number B ^ |l^| 

F(B) = ^ 2 B B 2 /V, (87) 

with us B = 13/6T 2 and V the volume. Due to the relation between a change in the Chern- 
Simons number to a change in the baryon number, the effective potential generates a force 
on the Chern-Simons number. As a result a non-zero baryon number will be "pushed back" 
to B = by the potential. This is the wash out of baryon number. 
In terms of O the effective potential is given by 

Kff(6) = ^(gvfcole 2 , (88) 

where the factor (gv) 3 is the inverse volume of a sphaleron. 

Now consider the effect of a first-order electroweak phase-transition. Besides the com- 
plicated dynamics of bubble nucleation and moving bubble walls, the quarks will acquire 
a mass through the Higgs mechanism. This produces a small change in the back-reaction 
of the quarks on sphaleron transitions. In the description of this back- reaction in terms of 
an effective potential, this means that the potential (|8^) is changed by mass corrections: 
uj 2 b —s- ll>b(1 — Cimf /T 2 ), with rrii the masses of the different quark species and Cj di- 
mensionless coefficients. As argued before, a change in a (harmonic) potential will lead to 
a temporary non-zero value of the one-point function. Hence, the expectation value of the 
baryon and Chern-Simons number will acquire temporarily a non-zero value. 

Before we discuss how "temporarily" may become forever, we will turn to the question 
of how large we may expect the expectation value of the baryon number to grow. From the 
sample calculations in Appendix C, we learn that the asymmetry has typically a maximum 
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value of ~ (1 — u>%/u>i)5, with S a dimensionless measure of the amount of CP-violation 
and Cl>i and Cl>2 the two frequencies of the potential before and after the phase transition. 
The factor (1 — ui^/uif), that is a measure of the amount of the departure from equilibrium, 
is of the order of the mass correction to the potential (j87|). Therefore, it is dominated by 
the top quark mass: (1 — oj^/ojf) ~ mj/6T 2 , with m t the top quark mass. In the case of 
Chern-Simons number diffusion we identify S with the relative asymmetry in the velocity 
at the sphaleron configuration Av rc i, S ph- Hence, when we consider the motion of O in this 
potential, we estimate for the maximum after the electroweak phase-transition 

9 

(e) max ~ ^|A Urel , sph |. (89) 

Imagine the universe made up of blocks of size (gv)~ 3 (the size of the sphaleron), with 
in each block a O-variable. At the first-order electroweak phase-transition bubbles fill out 
the universe and in each block the harmonic potential for changes. The time that the 
bubbles need to fill out the universe is much shorter than the time for a single sphaleron 
transition (g^) 3 !^^- Therefore, as far as the motion of O is concerned, the change in the 
potential occurs effectively at the same time in the whole universe. This implies that the 
expectation values of 0's in different regions in the universe move simultaneously towards 
their maximum value. Then we get for the maximum baryon-number density, n B = (B) /V , 



n B 



4^(6) max , (90) 



with the photon density n 7 = 0.24 T 3 . In this estimate we used that 0's in different blocks 
of space move independently, similar as in the case for pure diffusion in going from ( p7[ ) 
to ( |48| ) . This is only valid when the back reaction of the baryons on sphaleron transitions 
may be considered as arising from the potential ( |87| ) and non-localities may be neglected 
on the length scale {gv)~ x . This is what is usually done in the literature, see e.g. Q. If 
this assumption is incorrect extra suppression factors may be expected. 

When the generated baryon-number is not washed out, the asymmetry ( p0| ) results, at 
the present time, in an asymmetry 



n 7 



7#!p + 6#!p|10 4 f — — — 



CP f ""CP J 

max. now 



(91) 



where we used g = 0.6, m t = 170 GeV, and v ~ T ~ 100 GeV. This rough estimate for 
the asymmetry may indicate that in the process considered here, sufficient baryons may 
have been generated, provided that the CP-violation is strong enough and occurs at a not- 
too-high energy-scale. It may be noted that the mass M, that gives the energy scale of new 
physics, may not be as large as in scenario's for baryon generation based on CP-violation 
through CP-odd dimension-six operators, since (|9l]) is suppressed by four powers of M 
(instead of two) . Other constraints for the scenario based on asymmetric diffusion may be 
more important and are discussed in the next section. 

7.4 Constraints 

The immediate question is whether a once created baryon asymmetry survive to present 
times. The decay of the baryon-number density, ns, for small densities is given by H 



dt 

The formal solution is 



dn/3 o 

-9l3oj B T sph n B (92) 



n B {t) = n B (t pt )e- R{t \ (93) 
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with 



R(t) = 9 / dt'pu 2 B r sph {t'), (94) 



with t pt the time of the electroweak phase-transition and r sp h(i) the time-dependent 
sphaleron rate (its time-dependence enters through the changing temperature and Higgs 
expectation value as the universe cools down). The requirement that a generated baryon 
asymmetry survives to present times is 

R = R(oo)<l. (95) 

This may be translated into a bound on the phase transition strength 

v pt > v CI ~ 100 GeV, (96) 

with Vpt the Higgs expectation value just after the first-order electroweak phase transition. 

In |2(| one of the authors believed that the bound to prevent the wash out of baryon 
number could be much weaker. This was based on the realization that that the linear 
response result is not directly applicable to the situation where a non-zero expectation 
value of a coordinate is generated by asymmetric diffusion in a symmetric potential. (Since 
according to linear res ponse this expectation value can never become non-zero, which as 
was seen in section 7.1, is not the case.) However the time scale for equilibration given by 
linear response ( ^3| ) is still correct. Therefore, the bound ( p5| ) still applies. 

In the scenario discussed there is another bound, since the asymmetry itself has to be 
generated after the phase transition. Therefore, to generate a large asymmetry sufficient 
sphaleron transitions must occur after the phase transition. This requires 

R>1. (97) 

Hence, to generate sufficient baryons without a large wash out in the scenario that we 
discuss here, the number of sphaleron transitions after the phase transition has to lie in a 
small range around 1. We expect typically something like 

e" 1 < R < e. (98) 

This requirement for the number of sphaleron transitions after the electroweak phase- 
transition may be translated into a requirement on the strength of the phase transition, 
that may be indicated by the Higgs expectation value just after the transition, v pt . In 
appendix D we show that from ( p8|) it follows that this Higgs expectation value lies in the 
range 

/3pt-B S ph) pt 01 \ /3 p t-B S ph) ' ^ ^ 

with the inverse inverse temperature at the phase transition /3 pt ~ (100 GeV) -1 . Since, 
/Jpt-E-sph ~ 45 the range of allowed phase transition strengths is quite small. Hence, for the 
above scenario for the generation of matter to work, the strength of the electroweak phase 
transition has to satisfy strict bounds (^). For a given particle model these bounds may 
be translated into bounds for the Higgs mass, as has already been done for the bound ( |95| ) 
or ( |96| ) for the minimal supersymmetric standard model , p2| , p3| . 

7.5 Discussion 

Perhaps it is useful to consider the difference between the effect of the dimension-six oper- 
ator (j^4>FF (the lowest-dimensional CP-odd operator) and the dimension-eight operators 
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(D<fi)'(D(f>)FF and (FF)(FF) in the context of electroweak baryogenesis. As we have ar- 
gued the dimension-six operator does not affect the diffusion of the Chern-Simons number 
in equilibrium. However in a non-equilibrium situation it may have an effect. In particular, 
when the Higgs expectation value changes in time the dimension-six operator introduces an 
effective force F e s ~ dt{4>'<f>) on the motion of the Chern-Simons number. This implies that 
in the case of a time-dependent Higgs expectation value the diffusion of the Chern-Simons 
number yields a non-zero expectation value of iVcs itself: ([iVcs(*) — ^Vcs(O)]) ~ r s ph-Feffi- 
This in contrast to the asymmetry generated by the dimension-eight operators in equilib- 
rium that affect only expectation values of third and higher powers of Ncs{t) — iVcs(O). 

It is well-known that the dimension-six operator may play a role in the generation of 
the baryon asymmetry. For instance, during the electroweak phase-transition when the 
Higgs expectation value grows the effective force can be used to push the Chern-Simons 
number and baryon number to some positive value. The problem is that during the growth 
of the Higgs expectation value the sphaleron rate decreases. Therefore it is difficult to 
generate a sufficient baryon asymmetry in the short time that the Higgs expectation grows. 
A way to circumvent this problem has been given in |fj"5[ . fn that paper the situation 
was considered that the reheating temperature after inflation is below the temperature 
at which the electroweak phase-transition takes place. This means that there never was 
an electroweak phase-transition. The Higgs expectation value grows during the reheating 
of the universe. At this time the CP-odd dimension-six operator generates the earlier 
mentioned effective force. Further it has been argued in p5| that during the period of 
reheating the exponential suppression of the sphaleron rate is absent. Therefore sufficient 
sphaleron transitions can take place and a baryon asymmetry may be generated. Since the 
sphaleron rate becomes exponentially suppressed when the plasma thermalizes, a generated 
baryon asymmetry may be preserved. 

When the asymmetric diffusion of the Chern-Simons number from CP-odd dimension- 
eight operators plays a role in the generation of a baryon asymmetry, as we have argued is 
possible, then the time-scales for the generation and the wash out of a baryon asymmetry 
are similar. Therefore there is no need to try to avoid the exponential suppression of the 
sphaleron rate. The price to pay is for these similar time scales is that the strength of the 
phase transition needs to be finely tuned in order to ensure sufficient sphaleron transitions 
afte r the phase-transition and avoid a subsequent wash out of baryon number (see section 



Of course, in section |7| we have neglected various important dynamical aspects of the 
problem. Such as the motion of the bubble walls of a first-order phase transition, the motion 
of the baryons and dynamical aspects of their back-reaction on the motion of the Chern- 
Simons number, the growth in energy of sphaleron configuration as a bubble wall passes a 
certain region in space (that is as the Higgs expectation value increases) , the non-Brownian 
beginning of the motion of the Chern-Simons number (24|] etc. These neglected aspects 
may well modify the here presented estimate for the final baryon asymmetry. It is even 
possible that, being non-equilibrium processes, one of these will provide, in combination 
the included CP-violation, a different mechanism for baryon number generation. If so, the 
dimension-eight operators in (^) may prove more effective than the dim-six operator in ([}]) 
in providing the necessary CP-violation when the Higgs expectation value is not rapidly 
changing. Although all these non-equilibrium phenomena could yield a non-zero baryon 
number in combination with CP-violation, we believe that the mechanism presented in 
section (^) is the most natural way to generate an asymmetry on the basis of the CP-odd 
dimension-eight operators. 
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8 Summary 



We have argued that the inclusion of CP-odd operators, especially the dimension-eight 
operators in (Q), in an effective action will result in an asymmetric diffusion of the Chern- 
Simons number (in the absence of fermions). That is, in thermal equilibrium the correlation 
function ([Ncs(t) — Acs(O)] 3 ) becomes nonzero. In section @, we noted that, on general 
grounds, the expectation value of the third power of a coordinate either grows linearly in 
time or goes to a constant in the long time limit. We have argued that in the broken 
phase ([Ncs(t) — -/Vcg(O)] 3 ) grows linearly in time. A more detailed conjecture was given 
in (|48|). Unfortunately, the study presented here did not allow us to determine the sign of 
the expectation value (the direction of the growth) . 

Further we noted that asymmetric diffusion may lead to a non-zero expectation value 
of the coordinate itself in a non-equilibrium situation. For instance when the potential is 
changed (even if before, during, and after the change the potential is symmetric). This 
may have implications for the generation of the baryon asymmetry in the early universe. 
Although there are severe constraints on the strength of the electroweak phase-transition, 
we found that the observed baryon asymmetry could be generated by asymmetric diffusion 
in combination with the electroweak phase-transition. 
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A Relation between the Chern-Simons number and G 

The change in the Chern-Simons number is given by 

N CS (t) - A CS (0) = jL^dt' J d 3 xF^ a . (100) 

Using the parametrization (^|) of the non-contractable loop in configuration space found by 
Manton |l2), equation ( |l00|) may be rewritten in terms of O: 

Acs(t) - A cs (0) = - / *'esin 2 e / dr(d r f)f(l - f). (101) 

71" JO JO 

The time and spatial integrations may be factorized. The spatial integration is completely 
determined by the boundary values of the function / (it does therefore not depend on the 
Ansatz that is made for /, given the boundary conditions) 

j™ dr(d r f)f(l - /) = i (102) 

Hence, the relation between the change in Chern-Simons number and O is given by 

Ncstt) ~ A cs (0) = - f dt'Q sin 2 9. (103) 
TT Jo 

For instance when we change O from 0(0) = to O(t) = mr, crossing n barriers, the 
Chern-Simons changes as 

iVcs(*)-JVcs(0) = n, (104) 

as expected. 
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B Random walk model 



A random walk model is a convenient toy model to study diffusion. For asymmetric diffusion 
we consider a one-dimensional random walk model with probabilities 

P R = (1 + A)/2, P L = (1 - A)/2, (f 05) 

of moving right (i?) or left (L). The distance the particle moves in one time step differs 
also between the two different directions. When the particle moves to the right it travels a 
distance Ax R and when it moves to the left a distance Axl- 

Ax R = (I + A)- 1 , Ax L = (1- A)- 1 . (106) 

The relation 



PrAxr - P L Ax L = 0, (107) 

implies that the flux vanishes. 

Thus, the random walk model defined above contains the two important features of the 
motion of the Chern-Simons number established in section [| The flux vanishes (at every 
point), see (|39|), and the left and right velocities differ, see (ff5|), which translates in the 
different distances (106) travelled in one time step in the random walk model. 

The distribution function of the number of steps to the right, R and the number of 
steps to the left,L, after a total of N = L + R steps, reads 

P(L, R, N) = 6 L+R , N ^±(1 - A) L (l + A) R (f 08) 

When the particle starts at the initial time, N = 0, at position x — 0, x is given by 

x = RAxr — LAxl- (109) 



The relation L + R = N and ( |f 09| ) can be used to convert (108) in a distribution function 
of x. Instead, we will calculate expectation values. Let us start with (x). As mentioned 
before, the vanishing of the flux should imply that (x) is constant and, since the particle 
starts at x = Q, is zero. This may be checked for the random walk model. The expectation 
value of x is given by 

N N /VI 1 

L=0R=0 

x [R(l + A)' 1 - L{1 - A)- 1 ] . (110) 
We define R' = R—l and L' = L — 1 and write (pK 

N N-l 

" /.!/."! 2 v 



L=0 7?'=0 



N ~ 1 N ATI 1 

£ E ^ +R , N - ljmw {l - A) L (1 + A)* (111) 

L'=()R=0 



Since in the first line on the right hand side of (111) the term L = N in the sum over 
L does not contribute due to the delta function, and similarly the term R = N does not 
contribute in the second line, the sums in the first and second line cancel. Hence 

(x(N))=0, (112) 
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in accord with (|3G 

Next, consider the expectation value 

N N ATI 1 

<* 3 w> = EE*wzL^ (1 - A)i(1+A) * 

L=0 K=0 

x [i?(l + A)- 1 -L(l-A)- 1 ] 3 . (113) 
This expression may be evaluated in a similar manner as ( |110| ) . The result is 

(z 3 (N)) = \n [(1 + A)- 2 - (1 - A)- 2 ] . (114) 

Hence, for positive A, the expectation value (x ) is negative and grows linearly in time. 

This calculation supports our argument that when the asymmetry is independent of 
space and constant in time the expectation value (x 3 ) grows linearly in time. 

C Some solutions to the stochastic equation 

In this appendix we present some calculations using the Langevin equation ( |66| ) with a 
constant potential or a harmonic potential. 

We start with a calculation of (x) in the case V = 0. The Langevin equation ( J66] ) with 
V = reads 

x (1 + 3Sx + 6S 2 x 2 ) = - 7 x + f, (115) 

with 

(£(*)>« = 0, (116) 
(Z(tm')h = 2 7 TS(t-t% (117) 

and initial conditions x(0) = cc; n , x = Via- 

We solve the equations of motion in an expansion in 5: 

x(t) =x (t) +Sx!(t) + ..., (118) 

where xq and x\ satisfy the equations 

x a = -7A0 (H9) 

Xi = —jii — 3xoXo, (120) 

with initial conditions 

I (0) = X- m , ±0(0) = Win, (121) 

xx(0)=0, ±i(0)=0. (122) 

The solution for Xq reads 

x (t) = x in + — (1 - e-*) + f dt'G(t - t')t(t'), (123) 
7 Jo 

with the Green function G(t — t') = [1 — e _7( -* _t The solution for x\ is 

Xl (t) = -3 [ dt'G{t-t')x (t')x (t'). (124) 
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When the average over the stochastic force is performed, we obtain 

(x {t))t = a^ + ^i-e-r*), (125) 

7 

(*i(t)h = ^l-T]{\-e-^ + e-^). (126) 

We take thermal initial conditions: x- m — and v- m is thermally distributed. The average 
over initial momenta gives for the initial velocity and velocity squared 

(t*n}in = 0, (vf n )=T + 0(6). (127) 



Hence the average over initial conditions of ( 125 ) and ( 126 ) gives 

((xo(t))t)m = 0, (128) 
«*! (*)>£>>„ = 0{5). (129) 

Therefore the combined thermal and stochastic average of x vanishes, 

((x(t)) s ) ia = 0, (130) 

up to order 6. 

This is an explicit check that when the velocity is thermally distributed, and hence 
(x) = 0, that then (x(t)) = constant. 

The second quantity that we have calculated is (x 3 (t))^. We find that it goes to a 
constant (that depends on the initial distribution) in the long time limit. For the initial 
conditions x ln = and t>i n thermally distributed, we find 

«^W)«>in--^(^) 2 . (131) 

Hence, differently than we have argued for the diffusion of the Chern-Simons number, in 
this case the expectation value of the third power does not grow linearly in time. Since there 
are no high barriers, we don't expect the velocity to thermally distributed, independent 
of the position in space. Therefore the asymmetry in the velocities does not need to be 
independent of the position in space. In this way the argument for the linear growth of 
the expectation of x 3 (i) does not hold in this case. Hence, the fact that (x 3 (t)} goes to 
a constant when V = 0, does not contradict the expected linear growth of ([Ncs(t) — 
-^cs(0)] 3 ) m the long time limit. It does raise the following question. When we add a 
periodic potential VbSin 2 (x), how does the long time behavior of (x 3 (t)) depend on Vb? 
We know that for Vb = it goes to a constant, and we have argued that for Vb/T » 1 it 
grows linearly in time. So, at what value of Vb does the behavior change? 

Next, we consider the time-evolution in a harmonic potential V = \lo\x 1 . As initial 
distribution we use for the position exp — P^ujfxf n and v m is thermally distributed. This 
situation may be viewed as follows. The system is in thermal equilibrium before the initial 
time, t = 0, with a potential V — \uj\x 1 . At the initial time the harmonic potential changes 
^io\x 2 — ► \uj\x 2 . We calculate the following time evolution of (x). In the case without 
damping, 7 = 0, We get to first order in 6 

(x(t)) = —(l - ^l) sin(w 2 i)[l - cos(W)l. (132) 

For the case 7 >> u>, we present only the contribution of the slowest decaying mode 

( X (i)) = (1 _ 4) e~^h)t + 0{e -^l h )ty (133) 

27 V cjj / 



We learn from (132) and (133) that the expectation value (x) becomes non-zero after the 
potential is changed, as we have argued in section [7]. And further that the magnitude of 
(x) is proportional to (1 — LU^/ojf) and 5. 
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D Bound on the strength of the phase transition 

In this appendix we show how the bound ( p8[ ) 

er 1 <R<e, (134) 

with 

poo 

R=9 dt'/3uj 2 B r{t'), (135) 



may be translated into a bound on the Higgs expectation immediately after the electroweak 
phase transition, v pt . 

In order to simplify the following calculations we take k in (R) to be constant, then 

9/3c4r sph = —KTe-^-K (136) 

As a further simplifying approximation, we will assume that the Higgs expectation value 
stays constant after the phase transition and that the only time dependence enters through 
the temperature as [^H| 

T 2 = o.03^, (137) 

where M p \ = 1.2 10 19 GeV is the Planck mass. 

When the Higgs expectat ion value, v, is time-independent so is the sphaleron energy 
-E sp h, and the integration in (135) can easily be performed. The result is 

= ldj^L e -A**V* (138) 

with K — 1.08 k and (3 p t the inverse temperature at the electroweak phase-transition. 
It is useful to determine the value v p t — v cl - for which 

R=l. (139) 

Inserting the result ( |13S| ) for R and taking the natural logarithm gives 

lnR= In K^- - (3 pt E sph = 0. (140) 

^sph 

From |)| we obtain the value for n = 1.3 10 4 . Inserting this value for n in K and then in 
( p0| ) yields 

/3 P t£sph = 45. (141) 

This gives for the phase transition strength 

w CT ^ T pt - 100 GeV. (142) 

In |ll[] the time- or temperature-dependence of the Higgs expectation value is taken into 
account in the derivation of v CT . We are more interested in the window of phase-transition 
strengths determined by (134). 

A small shift of v p t by an amount of 6v, v p t = v CI + 5v, saturates the lower bound in 

dH 

R(v CI +Sv) =e~ 1 . (143) 
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This equation gives for 5v 



Vr 



(1 - (3 pt E sph ) = -1. (144) 



Thus the bound (134) translates into the following bound for the phase-transition strength 

Vpt 

V cr U - R I ) <Vpt <V CI (l + — , J, (145) 

for (3 pt E sph » 1. 
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